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Abstract 

Several aspects of the analysis of the data obtained with the cryogenic 
heat-and-ionisation Ge detectors used by the EDELWEISS dark matter search 
experiment are presented. Their calibration, the determination of their en- 
ergy threshold, fiducial volume and nuclear recoil acceptance are detailed. 



1 Introduction 



The dominant (~ 90%) component of the mass budget of the Universe may 
consist in Weakly Interacting Massive Particles (WIMPs), which could be the 
Lightest Supersymmetric Particles (neutralinos in most models) 1 . WIMPs 
would be present at the galactic scale as a halo of mass typically ten times 
larger than the visible part of the galaxy. The EDELWEISS collaboration has 
developped heat-and-ionisation Ge detectors j2] to measure recoils induced 
by elastic scattering of galactic WIMPs on a target nucleus. Constraints on 
the spin-independent WIMP-nucleon cross-section in the framework of the 
Minimal SuperSymmetric Model (MSSM) have been derived from the nuclear 
recoil rate measured with the EDELWEISS detectors 0, 111- We present in 
this paper the experimental details of these measurements. In Section we 
describe briefly the experimental setup and the method of detection of an 
energy deposit in the target. We then show the calibration procedure for the 
heat and ionisation signals (Section |3J), the trigger threshold determination 
(Section H|) and the tagging of the nuclear recoils (Section EJ). We finally 
present an original method to determine the fiducial volume of the detectors 
(Section EJ. 

2 The EDELWEISS detectors 
2.1 Experimental setup 

The experimental setup of the EDELWEISS-I experiment is described in [Hj 
and IQ. We simply recall that up to three detectors can be housed in a low 
background dilution cryostat working at a regulated temperature (27 mK 
in |3] and 17 mK in jl]). The EDELWEISS detectors are made of a ger- 
manium absorber (target for the incident particles) equiped with a thermal 
sensor and with metallic electrodes for charge collection. The simultaneous 
measurement of both phonons and charges created by a single interaction is 
therefore possible. 

The main characteristics of the detectors studied in this article are given in 
Table For all of these detectors, the absorber is a ~320 g Ge cylindrical 
crystal (~70 mm diameter and 20 mm thickness). Their edges have been 
beveled at an angle of 45° (Fig. Q). The electrodes for ionisation measure- 
ment are made of 100 nm Al layers sputtered on the surfaces after etching. 
The top electrode is divided in a central part and a guard ring, electrically 
decoupled for radial localization of the charge deposition. The bottom elec- 

1 See e.g. £Q for a review. 
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trode is the common reference. For the GGA1 and GGA3 detectors (GSA1 
and GSA3), a 60 nm hydrogenated amorphous germanium (silicon) layer 
was deposited under the electrodes in order to reduce the charge collection 
problems associated with events where the energy is deposited close to the 
detector surface. It has indeed been shown that the probability that charge 
carriers be collected on the same-sign electrode during the diffusion phase 
which preceeds the charge collection (dead layer problem) is reduced for this 
type of detectors [5, 6]. 

The thermal sensor consists of a Neutron Transmutation Doped germanium 
crystal (NTD) 2 , close to the metal-insulator transition. It is glued on a sput- 
tered gold pad near the edge of the bottom Al electrode (Fig. [TJ. The 
resistance of the DC-polarized GeA16 sensor was chosen to be ~ 3 MQ for 
GeA16 (Trunning ~27 hiK), and ranged from 3 MQ to 6 MQ at 17 mK for 
the other detectors. Reliable electrical contacts and heat links have been 
achieved by the ultrasonic bonding of gold wires (diameter 25 jttm) on gold 
pads. The thickness of these pads has been chosen to minimize the pro- 
duction of dislocations in the absorber caused by the bonding. A thermal 
analysis of the detectors will be published in Ref. j^j. 

2.2 Detection method 

The rise in temperature due to an energy deposit in the absorber gives rise to 
a variation AR of the thermal sensor resistance. When the sensor is polarized 
by a constant current I, AR then induces a voltage fluctuation AV across 
the resistor, which corresponds to the heat signal: 



The ionisation signal is obtained by collection of the electron-hole pairs cre- 
ated by the interaction in the germanium crystal polarized through a bias 
voltage applied to the electrodes. A low bias voltage 3 is required to limit the 
heating of the cristal due to the drift of the charge carriers, known as the 
Neganov-Luke effect [Zj. 

The energy Er deposited by a particle interacting in the detector can be 
determined by subtracting the Neganov-Luke effect from the heat signal : 



2 The NTD thermal sensors have been produced by Torre and Mangin for GeA16 and 
by Haller-Beeman Associates for the other detectors. 

3 During the data takings, the bias voltage applied to the top electrode varied from 
± 3 V to ± 9 V depending on the detector. 



AV — AR x I 



(1) 
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where V is the bias voltage and e 7 = 3 V the mean electron-hole pair cre- 
ation potential in germanium for 7-ray interactions (electron recoils). The 
variables Eh and Ei stand respectively for the heat and ionisation signal am- 
plitudes calibrated for 7-ray interactions following the procedure described 
in Section EHI 



We define the quenching variable Q as: 

C = |i (3) 

This variable is of particular interest in the case of WIMP search since nuclear 
and electronic recoils correspond to different ionisation efficencies. As Ei and 
Eh are calibrated using 7-rays, Q = 1 for electronic interactions by defini- 
tion. In the case of nuclear recoils (such as those that would be produced 
by WIMP interactions), this ratio is much lower: Q ~ 0.3. The simultane- 
ous measurement of heat and ionisation therefore provides an event-by-event 
identification of the type of recoils and thus gives an efficient method to re- 
ject the dominant 7-ray background. The precise definition of the rejection 
criteria is discussed in Sectional 

3 Calibration and resolution of heat and ion- 
isation signals 

3.1 Calibration of heat and ionisation channels 

The ionisation signal Ej is calibrated using a 57 Co source that can be in- 
serted in the liquid He bath of the cryostat to a distance of ~ 10 cm from 
the detectors, with only a ~ 0.5 cm thick copper shielding layer between the 
source and the detectors. The 122 and 136 keV peaks are clearly visible on 
the spectra (Fig. Efc), allowing a precise calibration of the ionisation signal. 
The linearity of the signal amplitude has been verified using the 46.52 keV 
line from 210 Pb (Fig. Eb) i n the detector environment and the 8.98 and 
10.37 keV lines from the decay of cosmic-ray induced long life isotopes 65 Zn 
and 68 Ge in the detector. The calibration factor is observed to be stable 
within a fraction of percent over periods of months. Because of the parasite 
capacitance between the centre and guard electrodes, a charge fully collected 
on an electrode also induces a signal on the other. This cross-talk of a few 
percents is purely linear and remains constant in time for a given detector. 
It can thus be easily corrected off-line (Fig. The ionisation signal Ej is 
defined as the sum of the guard ring and center electrode signal amplitudes 



3 



after correction of the cross-talk and calibration of the two channels. 
The heat signal amplitude Eh is periodically calibrated using the same 57 Co 
source. In contrast with ionisation, the heat signal appears to be very sensi- 
tive to long term drifts of the NTD temperature. It may for example vary by 
a few percent during several hours after transfers of cryogenic fluids. Between 
two 57 Co calibrations, the heat signal is therefore monitored on a continu- 
ous basis using the data from the low-background physics runs themselves 
by setting the average value of the Q ratio to 1 for electron recoils. The 
46.52 keV line from 210 Pb and the 8.98 and 10.37 keV lines associated with 
cosmogenesis activation of 65 Zn and 68 Ge in the detector (Fig. Et) are used 
to check the quality of the calibration of the heat signal. 
It should be stressed again at this point that the heat and ionisation signals 
are calibrated using 7-ray sources, which induce electron recoils. The Ej 
and Eh values thus correspond to the actual energy deposit for this type 
of interactions only, and are therefore expressed in keV electron equivalent 
(keV ee ). 

3.2 Resolution of heat and ionisation channels 

For each detector, the baseline resolutions of the heat and the two ionisation 
channels are regularly controlled through runs with an automatic random 
trigger. These runs show that the noises of the three channels are not corre- 
lated. The ionisation baseline resolution can therefore be written as : 

(°/) 2 = (^center) 2 + {^guardY ( 4 ) 

The 57 Co calibrations give a measurement of the resolutions for the ionisa- 
tion and heat signals at 122 keV. Typical values obtained for the detectors 
studied here are given in Table |2] 

We parametrize the heat and ionisation signals resolutions at a given electron- 
equivalent energy E as : 

a ItH (E) = ^K,) 2 + M) 2 (5) 

where the factors a/ and an are deduced from the resolution of the ionisation 
and heat signals at 122 keV. The resolutions of the 10.37 and 46.52 keV peaks 
observed in low-background physics runs fit well with the expressions ai t n{E) 
from Eq. (jSJ) (Fig. |2Jl). It can be noted that the resolutions at Ej ~10 keV ee 
-an energy below which most of the WIMPs signal is expected- is dominated 
by the baseline resolutions a® and a H - 

Finally, the recoil energy resolution can be computed from the heat and 
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ionisation signal resolutions using Eq. (J2J). The noises of both signals being 
uncorrelated, this resolution can be written as: 

^ = ^(i+£)W+g)W (6> 

In the case of GeA16, and for the bias voltage applied during the low- 
background physics run (V=6 V), the resolution values displayed in Table El 
lead to o Er ~ 8 keV FWHM around 30 keV. This value is reduced to 4 keV 
FWHM in the condition of the low-background physics run recorded with 
GGA1 (V=A V) g]. 

4 Threshold 

The ionisation and heat channel data are continuously digitized and filtered 
at a rate of 200 kHz and 2kHz, respectively. When a filtered ionisation value 
exceeds a fixed threshold value, data samples in all detectors are stored to 
disk. The trigger is defined by requiring a minimum threshold on the absolute 
value of any of the filtered ionisation channels. For each event, the list of all 
detectors having triggered is stored bit pattern. 

The ionisation threshold value, Ei^h is defined as the ionisation energy (in 
keV ee ) at which the trigger efficiency reaches 50%. It is the most important 
parameter governing the recoil energy dependence of the efficiency. Its value 
is measured using two different techniques: one is based on the Compton 
plateau observed with a 7-ray source, and the other on coincidence neutron 
data. 

In the first one, a 7-ray spectra is recorded using a source producing a 
important Compton plateau, such as 60 Co or 137 Cs. Monte Carlo simula- 
tions indicate that the shape of the plateau above 10 keV can be linearly 
extrapolated to lower energy. The efficiency as a function of Ej, e(Ej), is 
thus obtained by dividing the measured rate by the straight line extrapolated 
from the rate above 10 keV. The resulting e(Ei) data is fitted by a integral 
of a gaussian (erf), yielding the experimental value of E I>th . However, this 
method is limited by the large data sample necessary to obtain a significant 
number of events in the threshold region. 

The second technique was made possible by the simultaneous operation 
of three detectors with a 252 Cf neutron source (and thus could not be ap- 
plied to the GeA16 detector). Neutron scattering induces a large number of 
coincidence events where at least two detectors are hit. The upper pannel of 
Fig. |U shows the Ej distribution recorded in one detector with the condition 
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that any of the other two detector triggered (unfilled histogram). Despite 
that the detector under study is not requested in the trigger pattern, the 
peak at Ei=0 due to baseline noise is not overwhelmingly large, due to the 
importance of the coincident rate. An unbiased sample of events with Ej > 2 
keV is thus obtained. When in this sample it is further requested that the 
detector under study be present in the trigger pattern, the shaded histogram 
is obtained. The ratio of the two distributions shown in the lower pannel of 
Fig. |3] correspond to the efficiency e(Ej). This interpretation is valid in the 
region close to Ei >t h an d above because in that energy range the contribu- 
tion of the peak due to baseline events is negligible and because the slope of 
the unbiased distribution is reasonnably small compared to the experimen- 
tal resolution on Ej. Indeed, applying this method to a distribution N(Ej) 
proportionnal to exp(— E^/r) and smeared with an experimental r.m.s. res- 
olution a, this method would result in a shift of — cr 2 /r of the deduced value 
of Ei >t h relative to the true value. In the present case, where the range of 
exponential slopes and resolution are 3 < r < 8 keV and 1 < 2.35a < 2 keV, 
the shift should not exceed 0.2 keV. 

Both Compton and neutron coincidence techniques give consistent ioni- 
sation threshold measurements. The coincidence measurements are the most 
precise, as the neutron source has the advantage of yielding a maximum rate 
at the lowest energy, and in addition, the quenching of ionisation for nuclear 
recoils ensure that the stability of the measurement can be tested by im- 
posing a cut on the heat signal Ejj without affecting the ionisation signals 
with Ei above ~ Eh/2. The measured Ej^h values for the different ionisa- 
tion channels of the detectors under study are listed in the last column of 
Table El 

5 Nuclear recoil band 

Figure El shows a (Er, Q) distribution from the data recorded with a 252 Cf 
source emitting 7-rays and neutrons. Experimentally, the Q variable appears 
to follow a gaussian distribution at the ~ 2a level for both nuclear and 
electron recoils populations (Fig. We therefore parametrize the region of 
90% acceptance for the nuclear recoils by the following cut: 

\Q- < Qn > I < 1.65<x Qn (7) 

where < Q„ > and ag n are the average value and the standard deviation 
of the Q distribution for nuclear recoils, both variables being determined 
for each detector from 252 Cf calibration data under the same experimental 
conditions as the low-background physics runs. 
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5.1 Neutron line 



The neutron line is the average Q value for the nuclear recoils population. It 
is parametrized from 252 Cf calibration data by : 

<Q n >(E R ) = a(E R ) b (8) 

The a and b values resulting from the fit of the experimental data for each 
EDELWEISS detector are statistically consistent with the values a = 0.16 
and b = 0.18 quoted in |S]. The biases on the determination of < Q n > 
due to experimental calibration uncertainties, heat quenching effects , and 
multiple scatterings are globally taken into account with this measurement. 

5.2 Electron and nuclear recoils zones standard devi- 
ations 

The standard deviation of the electronic and nuclear recoil distributions, 
respectively noted <tq 7 and crg n , can be calculated with Eqs. (j2J) and (j3J) by 
propagation of the experimental values 07 and cr#: 

^oJEr) = ^( 1 + i <Qn> ) + + <Qn> 2 ^m 

In the case of 60 Co, 252 Cf calibrations and low-background physics runs, the 
experimental values of <7q i at high energy are significantly larger (up to 
~+30% at 122 keV) than those calculated from the resolutions given in 
Table El with Eqs. ©, © and © (Fig. Ek). A dependance of the heat 
signal amplitude on the position of the interaction provides an explanation 
for this discrepency. This hypothesis is consistent with the ~ 1% heat sig- 
nal amplitude difference observed between center and guard events in 57 Co 
calibrations. We therefore enlarge the an coefficient in Eq. (JHJ so that the 
analytic expression given in Eq. © for aQ j (Ej{) actually follows the exper- 
imental distribution for 60 Co, 252 Cf and low-background physics runs. We 
have checked that 90% of the experimental events then fall inside the elec- 
tron recoil zone defined in this way. 

Even after correcting the a# value, the nuclear recoils Q distribution of 252 Cf 
calibration data is broader at high energy than what is expected from Eq. 
(fTTH) (Fig. Ub)- Atomic scattering processes [TT], fluctuations in the number 
of charges created by a nuclear recoil [TJ] and multiple scattering (see Section 



7 



15.31) are in particular expected to give an intrinsic width to the Q distribu- 
tion for nuclear recoils and thus explain this behavior. The experimental 
<jQ n dependance on recoil energy is properly described when a constant C is 
quadratically added to the term associated with the experimental resolution. 
The equation (|TU|) is thus re-written as follows: 

* Qn (E I d = y/o Qn (E R ) + C* (11) 

Typical values of C ~0.040 are determined for each EDELWEISS detector 
by fitting the experimental OQ n points using Eq. pijl . With this definition, 
we have checked for each detector that 90% of nuclear recoils induced by 
252 Cf calibrations are inside the nuclear recoil zone defined in Eq. (|7|). 

5.3 Effect of multiple scattering 

The nuclear recoil zone is determined through neutron calibrations, for which 
the proportion of multiple interactions is around 40% between 20 and 200 keV. 
This is of particular importance because in contrast to neutrons, WIMPs are 
expected to interact only once in the detector, and the Q variable is in this 
case larger than when the same energy is deposited in multiple nuclear in- 
teractions, as can be deduced from Eq. (JH1). 

We therefore evaluated quantitatively the effect of multiple interactions using 
a GEANT [T0| simulation of 252 Cf calibrations of the EDELWEISS detectors. 
The Q variable has been calculated for the simulated nuclear events by as- 
sociating with Eq. (jHJ) an ionisation signal of amplitude e/ = 0.16 (e^) 1 ' 18 to 
an energy deposit cr in a single interaction, and summing each individual 
ej to obtain the total Ej energy for a given neutron. The effect of multiple 
interactions has then been evaluated with these simulated data by smearing 
the resulting Q distribution with the experimental resolution given in Eq. 
(jlO)) . and then by comparing the distributions obtained when selecting or 
not single interactions events (Fig. |HJ). 

Although multiple interactions tend to lower < Q n >, this effect remains 
weak, and the Q distribution associated with single interactions events is 
only slightly narrower and completely included in the wider band. The nu- 
clear recoils zone determined through 252 Cf calibrations has therefore been 
conservatively used for the low-bakground physics run analysis. 

5.4 Analysis energy range 

Equations © and (|T0|) predict that the discrimination between electronic 
and nuclear interactions is deteriorated at low energies (see also Fig. 
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Rejection of the 7-ray background at a given level therefore defines a lower 
bound for the analysis energy range. 

Secondly, the detection efficiency has to be as close to 100% as possible in 
the analysis window in order to insure a good quality for the data set. The 
trigger threshold is therefore another factor which has to be taken into ac- 
count for the definition of the analysis lower energy bound. For both 2000 [3] 
and 2002 [3] runs, the choice of the analysis lower bound has mainly been 
driven by this last factor. The threshold values of 5.7 and 3.5 keV ee for the 
ionisation signal indeed correspond respectively to recoil energies of 30 and 
20 keV for a 100% detection efficiency, and 90% efficiency when the nuclear 
recoil zone is taken into account. 

Extensive 7-rays and neutron calibrations are performed before the physics 
data taking is initiated in order to fix the lowest recoil energy value corre- 
sponding to acceptable levels of 7-ray background rejection and detection 
efficiency. This ensures that the lower limit of the analysis window is not 
influenced by the possible presence of events in the final data set. The defi- 
nition of the upper bound of the analysis window is described in jlj. 

6 Fiducial volume 

6.1 Modelisation of the collection process 

The segmentation of the upper charge collection electrode in a central part 
and a guard ring leads to the definition of a fiducial volume. This volume 
is shielded against a significant amount of the radioactivity of the detector 
environment by the peripherical volume, as shown in [Hj. To allow for the 
experimental resolution on the ionisation signals, the fiducial cut is defined 
as corresponding to a fraction of 3/4 of the charge collected on the center 
electrode. In order to give a robust and precise estimation of the detector 
volume associated with this fiducial cut, it is necessary to relate a given ratio 
of the two ionisation signal amplitudes to a given volume inside the detec- 
tor. This is not a straightforward process: first, for non-WIMP interactions, 
multiple interactions have to be taken into account, and furthermore, inter- 
actions between charges may play a crucial role in the collection process. In 
particular, the important proportion of events with a charge signal shared 
between the two channels observed in each detector for 60 Co calibrations (see 
e.g. Fig. 3) hints to the importance of these charge interactions processes. 
In order to test their influence on the determination of the fiducial volume, 
we choose to model the collection process with the simplified phenomenolog- 
ical description of charge collection given in Ref. [T3], associated with the 
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hypothesis of a plasma effect before charge drift. We will see that, even if 
some of our results cannot be explained in the framework of this very sim- 
plified model (Section 16. 2|) . it provides a good empirical tool to determine 
the fiducial volume and estimate systematic errors on its value (Section l6.3j) . 
The model used here assumes the distribution of the charges in a sphere with 
uniform density, extending to a maximal radius r b before the charge is fully 
collected. Charges are distributed among the two electrodes depending on 
the position of the interaction relative to the surface corresponding to the 
separation between drift lines going to the center and guard rings. Here, we 
assume for simplicity that this surface is parametrized by a cylinder of radius 
Rc (Fig. fTTj) . For an interaction at the radius R > Rc + r b in the crystal, 
the whole charge is fully associated with the guard ring. If R < R c — r b , 
then the charge has to be associated with the center electrode. Finally, if 
Rc + r b > R > Rc~ r b, then the charge is splitted among the two electrodes, 
with a relative proportion associated with the center electrode correspond- 
ing to the fraction of the sphere inside the cylinder of radius Rc- For given 
values of r b and Rc, the fiducial volume is determined in this model by the 
following expression of the fiducial radius: 

Rfid = Rc + 2 cos [j^j r b (12) 

A fraction of 1/4 of the total volume of a sphere of radius r b centered on 
Rfid is inside the cylinder of radius Rc- In the framework of our model, 
interactions inside the cylinder of radius Rf^ thus correspond to a charge 
collection equal or greater than 3/4 of the total charge. 



6.2 Validity and limits of the modelisation 

In order to test its ability to reproduce the distribution of charge amplitudes, 
ionisation signals are simulated in the framework of this simple model, us- 
ing the program GEANT ^Hj for 60 Co and 252 Cf calibrations, as described 
in Section 15. HI The parameters r b and Rc of the simulated data are then 
adjusted to match the experimental distribution of the Y variable on a given 
energy range, the Y variable being defined as the normalized difference of 
the ionisation signals: 

Y guard -^center (13) 

TP I IP \ / 

^ guard i" ^center 

The result of this optimisation is shown in Fig. El in the case of a 60 Co cal- 
ibration of the GeA16 detector under 6.3 V bias voltage. The shape of the 
simulated distribution closely follows that of the experimental data, while 
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a simulation using an alternative model (linear distribution of the charge, 
detailed in clearly exhibits a different pattern. 

We have also studied the evolution of the (r&, Rc) parameters as a func- 
tion of bias voltage for the GeA16 detector [THj. Rc should not depend on 
the value and sign of the bias voltage, since it is related to the static field 
distribution only, while should increase with decreasing bias voltage: as 
the field increases, the less time there is for diffusion processes. The values 
of the parameters r& and Rc determined for 60 Co and 252 Cf calibrations of 
the GeA16 detector versus the applied field are displayed in Fig. HH The 
and Rq values follow the expected behavior. Moreover, the mean measured 
value of Rc (< Rc >— 24.45 ± 0.05 mm) is statistically compatible with the 
value Reiectro — 24.4 mm expected from numerical calculations of the electric 
potential inside this detector (see Table HJ) . 

The very large r& values are a clear sign that macroscopic charge extension 
perpendicular to the drift direction occurs before the charge collection is 
completed. However the data does not support that this expansion is driven 
by the plasma effect invoqued in Ref. [12]: a charge cloud size of the order 
or above a millimeter is indeed not compatible with results of studies on the 
dead layer [5, 16]. Furthermore, Fig. ITUl shows that r&(— ) << rf,(+) and 
that the values of r& for 60 Co and 252 Cf calibrations do not differ significantly 
for a same bias voltage. These two experimental results are also in strong 
disagreement with the predictions derived from the hypothesis of a plasma 
effect: firstly, the observed asymmetry for r& values between positive and 
negative bias voltage does not find any explanation in the framework of the 
plasma model, and secondly, the plasma effect should be weaker in the case 
of 60 Co calibrations than for 252 Cf (and thus values much smaller), since 
7-rays induce much lower charge densities than neutrons. These are strong 
indications that the simple model presented here does not provide a proper 
description of the dynamics of the charge drift and collection. Charge re- 
pulsion during drift, not taken into account here, could for example play an 
important role in the collection process. A more detailed study, with dedi- 
cated detectors, has been initiated in the EDELWEISS collaboration in the 
aim of better understanding the collection process |16j . 

6.3 Measurement of the fiducial volume 

Our results clearly point out the limits of the modelisation presented in Sec- 
tion 16.11 Still, it has to be stressed that this model reproduces correctly the 
distribution of charges among the electrodes (Fig. which represents the 
net effect of the charge collection process. It is therefore sufficient to give a 
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precise determination of the fiducial volume and evaluate possible systematic 
errors, before a better, physically motivated model replaces it. 
We have calculated Rfid with Eq. |T2|) and the and Rc values deter- 
mined from 252 Cf calibrations under the same bias voltage as that of the low- 
background physics run for each detector. For all detectors except GeAllO, 
the Rc values are compatible with those expected from the geometry of the 
electrodes and from numerical simulations of the electric potential inside the 
detectors. The values determined for Rfid for the detectors are summarized 
in Table El The systematic error associated with the uncertainty on the ex- 
act mechanism producing the charge expansion is evaluated by taking the 
difference between the fiducial volume value deduced using the linear model 
and the one presented in Section 16.11 Despite the poor description of the 
charge distribution by the linear model (see Fig. EJ), this difference is only 
1%. The variation of the energy range used to determine the values of and 
Rc through comparison of the experimental and simulated Y distributions 
proved to be a minor contribution to this systematic error. 

An alternative evaluation of the fiducial volume is the fraction of cosmic 
activation events at 8.98 and 10.37 keV (see Fig. EJ) selected by the fidu- 
cial volume cut. Such events are expected to be evenly spread inside the 
detector, and are observed at rates varying between 3 to 15 events per detec- 
tor per day. In the few days following a neutron calibration, the 10.37 keV 
rate is also enhanced due to n Ge activation (Ty 2 = 2.7 d), a population 
that is also expected to be evenly spread inside the detector. The measured 
fractions, directly interpreted as Vfia values, are listed in the last column of 
Table El They are compatible within statistics with the values derived from 
the neutron calibration data and the collection process modelisation. The 
cosmic activation data is however less precise due to statistics, but this mea- 
surement is a good cross-check for the determination of the fiducial volume, 
and validates the use of the model presented in Section l6~Tl to determine the 
fiducial volume. 

7 Conclusion 

We have described in the present work the calibration aspects of the data 
analysis in the EDELWEISS experiment. In particular, the nuclear recoil 
zone and fiducial volume have been estimated using several methods, allow- 
ing to define a conservative value of these important parameters. A simple 
parametrization allows us to reproduce accurately the distribution of the 
charges between the centre and guard electrodes associated with 60 Co and 
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252 Cf calibrations, making possible the systematic studies necessary to es- 
tablish the robustness of the determination of the fiducial volume of the 
detectors. 
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Label 


Mass 
(g) 


D 

(mm) 


Vol. NTD 

(mm 3 ) 


Amorphous 
layer 


T 
(mK) 


GeA16 


321.62 


24.4 


4.0 


none 


27 


GeA19 


325.43 


24.0 


5.6 


none 


17 


GeAllO 


323.91 


24.0 


5.6 


none 


17 


GGA1 


318.50 


24.0 


1.64 


Ge 


17 


GGA3 


324.40 


24.0 


5.6 


Gc 


17 


GSA1 


313.68 


24.0 


5.6 


Si 


17 


GSA3 


297.03 


24.0 


5.6 


Si 


17 



Table 1: Main parameters for the EDELWEISS detectors studied in this 
article. "R e iectro " refers to the radius value of the cylindrical volume associ- 
ated with charge collection on the center electrode. These values are calculated 
through electrostatic simulation of the detector, taking into account the actual 
electrodes geometry. The existence of an amorphous Ge or Si layer under the 
electrodes is also mentionned. "T running " is the value of the regulated cryostat 
temperature while running. 
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FWHM @ keV 


FWHM 


@ 122 keV 


Trigger Threshold 




Center 


Guard 


Heat 


Ion. 


Heat 


Center 


Guard 


Detector 


(keV ee ) 


(keV ee ) 


(keV ee ) 


(keV ee ) 


(keV ee ) 


(keV ee ) 


(keV ee ) 


GeA16 


2.0 


1.4 


2.2 


2.8 


3.5 


6.0 


4.0 


GeA19 


1.2 


1.4 


0.5 


2.6 


3.3 


4.3 


4.9 


GeAllO 


1.1 


1.3 


0.4 


3.0 


3.5 


3.3 


4.3 


GGA1 


1.3 


1.3 


1.3 


2.8 


3.5 


3.5 


3.5 


GGA3 


1.3 


1.5 


0.4 


3.1 


2.7 


2.9 


3.9 


GSA1 


1.2 


1.4 


0.6 


3.1 


2.8 


3.5 


3.4 


GSA3 


1.1 


1.3 


1.4 


3.3 


3.3 


3.0 


3.4 



Table 2: Typical values obtained in keV ee for the full width half maximum 
resolution for heat and ionisation signals at and 122 keV for the detectors 
studied in this article. The precision on these measurements are ±0.1 keV 
at keV and ~ ±0.2 keV at 122 keV. Also given here are the threshold 
values for the center and guard channels. The precision is ±0.1 keV for both 
channels, except for GeAW where it is ±0.5 keV. 



Detector 


Bias 


n 


Rc 


Rfid 


Vfid 


Activation 




(V) 


(mm) 


(mm) 


(mm) 


(%) 


v fid (%) 


GeA16 


+6.34 


4.3 + 0.2 


24.5 + 0.3 


23.0 + 0.3 


54.6 + 1.4 


50 + 3 


GeA19 


+2.00 


6.1 + 0.2 


23.5 + 0.4 


21.4 + 0.4 


47.4 + 2.0 


53 + 4 


GeAllO 


-3.00 


2.3 + 0.2 


21.9 + 0.4 


21.1 + 0.4 


46.0 + 1.7 


50 + 4 


GGA1 


-4.00 


1.6 + 0.1 


24.6 + 0.2 


24.1 + 0.2 


60.1 ± 1.1 


57 + 3 


GGA3 


-4.00 


1.4 + 0.1 


24.1 + 0.2 


23.6 + 0.2 


57.7 + 0.7 


60 + 5 


GSA1 


-4.00 


1.5 + 0.1 


24.2 + 0.2 


23.7 + 0.2 


58.3 + 0.8 


61 + 4 


GSA3 


-4.00 


1.5 + 0.1 


23.9 + 0.2 


23.3 + 0.2 


56.2 + 0.8 


61 + 5 



Table 3: Values of various parameters for the EDELWEISS bolometers 
determined from 252 Cf calibrations under the given bias voltage. The error 
bars correspond to statistical errors. The systematic error on Vfid is 
"Activation" refers to the fraction of 8.98 and 10.34 keV events recorded with 
the fiducial volume cuts (E center > 3E guard ). 
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Figure 1: Top pannel: EDELWEISS GGA1 detector ($ =70 mm). Bottom 
pannel: close-up on the NTD thermal sensor glued on its golden pad on the 
beveled part of the crystal. 
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Ionization - centre (keV) Heat signal (keV) 

Figure 2: Figs, a), b), c): spectra obtained with the GeAUO and GGAl 
detectors in various energy regions during physics background runs ( a and b ) 
and 57 Co calibration runs (c). For the first two spectra, the lines associated 
with 65 Zn (8.98 keV) and 68 Ge decays (10.37 keV) (a) and 210 Pb contami- 
nation (46.52 keV)(b) are clearly visible. In Fig. d), the baseline and peak 
resolutions of the GeAUO detector heat channel for the 10.37, 46.52 and 122.1 
lines are fitted by the expression given in Eq. J3J). 



18 




Figure 3: Distribution of the guard ring versus the center electrod signals for 
a 60 Co calibration of the GeAW detector at a bias voltage of +6.34 V before 
(a) and after (b) linear correction of the cross-talk between the two channels. 
On Fig. (b), the events along the horizontal (vertical) axis correspond to 
center (guard) events, for which the charge is fully collected on the center 
electrode (guard ring), and the events between the two axis correspond to 
shared events. Shared events represent a proportion of ~ 50% of the total 
number of events for this calibration. 
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Figure 4: Top pannel: plot of the data recorded in the GGA1 detector dur- 
ing neutron calibration for any other detector (GeAW or GeAUO) triggering 
(line), and with the additional condition that GGA1 also triggers (hatched 
area). Bottom pannel: experimental efficiency curve for the GGA1 detector 
corresponding to the ratio of the two distributions from the top pannel. The 
100% efficiency is reached at 5.5 keV ee energy, corresponding to 20 keV recoil 
energy. 
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E R (keV) 



Figure 5: Projection in the (Er, Q) plane of the events recorded in the 
GGAl detector during a 252 Cf calibration. The thick lines represent the 90% 
nuclear and electronic recoils zone (=tl.645cr around < Q„ > and < Q 7 > 
respectively). The dotted line corresponds to the ionisation threshold curve 
(Ei=3.5 keV ee in this case). The dashed lines show where events associated 
with the inelastic scattering of neutrons on 73 Ge (13.26 and 68.75 keV excited 
levels) are expected in this plane. 
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Figure 6: Top pannel: spectrum of the Q variable in the 40-200 keV re- 
coil energy range for events recorded in the GGA1 detector during a 60 Co 
calibration (electron recoils). No events are seen for Q < 0.6. This shows 
the excellent quality of the charge collection for this detector. This test is 
performed for every detector before a low-background physics run is started. 
Bottom pannel: Spectrum of the Q variable in the 40-60 keV recoil energy 
range for events recorded in the GGA1 detector during a 252 Cf calibration 
(nuclear and electronic recoils). As the fit shows, the nuclear and electron 
recoils populations follow gaussian distributions down to the 2a level. 



22 




Figure 7: Experimental values for <tq (a) and OQ n (b) for a 252 Cf calibra- 
tion of the GGA3 detector. Also shown are the computed laws for these two 
variables from Eqs. (0|) and 177]) . before (dotted line) and after (solid line) 
correction of the an factor for (a), and for C=0 (dotted) and 0=0.0^0 (solid) 
for (b). The large value of ctq 7 in (a) in the 60-80 and 80-100 keV recoil en- 
ergy bins (bright color) is due to the inelastic scattering events (Fig. [3J) ; 
while the low value of OQ n in (b) in the 10-20 keV energy range (bright color) 
is related to the ionisation threshold, which artificially narrows the nuclear 
recoils distribution for low energies. The vertical dashed lines correspond to 
the bins limits. 
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Figure 8: Normalized spectra of the Q variable for a GEANT simulation 
of a 252 Cf calibration, selecting (hatched area) or not (continuous line) single 
interactions, in the 40-60 keV (top pannel) and 90-200 keV (bottom pannel) 
energy ranges. The distribution corresponding to all interactions is slightly 
shifted down (~ 0.015 units between 20 and 200 keV) in regards to the single 
interactions distribution, and is only slightly broader. 
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Figure 9: Distribution of the Y = Eauard ^ center variable for events of 

& guar di-fc' center 

60 Co calibration under +6,3 V bias voltage with 30 keV< Ej <1400 keV 
(cross). The two other spectra correspond to the simulated distributions ob- 
tained for modelisations in the case of a charge reparted in a sphere (dashed, 
rb=4-2 mm, Rq=24-5 mm) and charge reparted linearly (dotted, for param- 
eters giving the same fiducial volume). 
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Figure 10: Optimized values for rj, (top pannel) and Rc (lower pannel) for 
60 Co (circles) and 252 Cf (triangles) calibrations of the GeAW detector versus 
applied field values. Positive (negative) bias voltages correspond to the empty 
(filled) symbols. The r^ distribution is fitted by a power law: r^ = aE~ b . 
The fit gives a+ = 5.7 ± 0.1, a- = 1.8 ± 0.1, b + = -0.24 ± 0.02 and b~ = 
-0.28 ±0.09. 
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Figure 11: Representation of the GGAl detector in the (R,Z) plane. The 
thick line corresponds to Rfid = 24.1 mm. Also shown is the line Rc = 
24.6 mm delimiting the volumes associated with collection on the center and 
guard electrodes (dotted) and a circle of radius r^=1.6 mm centered on Rq 
(dotted). These values correspond to the optimum for a 252 Cf calibration of 
the GGAl detector under a bias voltage of -4-00 V. 
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